Possible CP-violation effects in core-collapse Supernovae 
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We study CP-violation effects when neutrinos are present in dense matter, such as outside the 
proto-neutron star formed in a core-collapse supernova. Using general arguments based on the 
Standard Model, we confirm that there are no CP-violating effects at the tree level on the electron 
neutrino and anti-neutrino fluxes in a core-collapse supernova. On the other hand significant effects 
can be obtained for muon and tau neutrinos even at the tree level. We show that CP violating effects 
can be present in the supernova electron (anti)neutrino fluxes as well, if muon and tau neutrinos 
have different fluxes at the neutrinosphere. Such differences could arise due to physics beyond the 
Standard Model, such as the presence of flavor-changing interactions. 
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I. INTRODUCTION 



■ Recent results from solar, atmospheric and reactor experiments have significantly improved our knowledge of the 
' neutrino mass differences and of two of the mixing angles. If the remaining mixing angle, #13, is relatively large there 

is a possibility that violation of CP symmetry may be observable in the neutrino sector. Currently planned and future 
experiments will have improved sensitivities to the value of this angle (see e.g. 0,0, ML, Effects of CP-violation in 
, accelerator neutrino oscillation experiments have been extensively investigated [ElElEHIEil- The discovery of a non 

■ zero Dirac delta phase might help our understanding of the observed matter-antimatter asymmetry of the universe 
CN \ HE EH E3) EH • Besides studies on terrestrial experiments with man-made sources, a few recent works have addressed 

CP-violation with neutrinos from astrophysical sources (see e.g. [3, EH)- The purpose of the present paper is to 
explore possible effects coming from the CP-violating phase in dense matter, such as that encountered in core-collapse 
supernovae. 

Core-collapse supernovae occur following the stages of nuclear burning during stellar evolution after an iron core 
is formed. The iron cores formed during the evolution of massive stars are supported by the electron degeneracy 
pressure and hence are unstable against a collapse during which most of the matter is neutronized. Once the density 
exceeds the nuclear density this collapse is halted. Rebounding pressure waves break out into a shock wave near 
the sonic point where the density reaches the nuclear density. Evolution of this shock wave and whether it produces 
an explosion is a point of current investigations. However, it is observed that the newly-formed hot proto-neutron 
star cools by neutrino emission. Essentially the entire gravitational binding energy of eight or more solar mass star 
is radiated away in neutrinos. Although the initial collapse is a very orderly (i.e. low entropy) process, du ring the 
cooling stage at later times a neutrino-driven wind heats the neutron-rich material to high entropies [IE EE Ell- 
Neutrino interactions play a very important role in the evolution of core-collapse supernovae and in determining 
the element abundance Neutrino heating is a possible mechanism for reheating the stalled shock |2fj|. A good 
fraction of the heavier nuclei were formed in the rapid neutron capture (r-process) nucleosynthesis scenario [2l| . Core- 
collapse supernovae are one of the possible sites for the r-process nucleosynthesis. A key quantity for determining the 
r-process yields is the neutron-to-seed nucleus ratio, determined by the neutron-to-proton ratio, which is controlled 
by the neutrino fluxes. In addition, recent work indicates that neutrino-neutrino interactions plays a potentially very 
significant role in supernovae [H |H [H, HE HE H3] . 

In this paper we study CP violation aspects in the core-collapse supernova environment. We first analyze analytically 
and in general terms, how the neutrino propagation equations and the evolution operator are modified in matter, in 
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presence of a non-zero Dirac delta phase. We obtain a general formula which is valid for any matter density profile 1 . 
In particular we demonstrate that, as in vacuum, the electron (anti)neutrino fluxes are independent of the phase 5, if 
mu and tau neutrinos have the same fluxes at the neutrinosphere in the supernova 2 . On the other hand the electron 
(anti)neutrino fluxes will dep end on S, if mu and tau neutrinos have different fluxes at the neutrinosphere, at variance 
with what was found in [14j |. We present numerical calculations on possible CP violation effects on the mu and tau 
neutrino fluxes as well as on the electron (anti-)neutrino flux and the electron fraction. The latter can only appear 
if physics beyond the Standard Model, such as flavor changing interactions, induces differences on the mu and tau 
neutrino initial total luminosities and/or temperatures. Finally we calculate the effects from the CP violating phase 
on the number of events in an observatory on earth. 

The plan of this paper is as follows. In Section II we present the general formalism to describe the neutrino 
evolution in presence of the S phase. The formulas concerning neutrino fluxes and the electron fraction in the 
supernova environment are recalled in Section III. Numerical results on these quantities are presented in Section IV. 
Conclusions are made in Section V. 



II. NEUTRINO MIXING IN THE PRESENCE OF CP- VIOLATING PHASES 



A. Neutrino mixing in ordinary matter in presence of CP-violating phases 



The neutrino mixing matrix is U a i where a denotes the flavor index and i denotes the mass index: 



(1) 



For three neutrinos we take 

U a i — 723T 13 T 12 = 




C 13 S 13 e~* s 

1 
Si 3 e* s C 13 




(2) 



where C13, etc. is the short-hand notation for cos #13, etc. and S is the CP-violating phase. The MSW equation is 
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for the charged-current and 



V c (x) = V2G F N e (x) 



V n {x) = --j=G F N n (x). 



(3) 



(4) 



(5) 



for the neutral current. Since V n only contributes an overall phase to the neutrino evolution we ignore it 3 . Following 
references (30| and [3l| we introduce the combinations 



cos6»23* M -sin6' 2 3l' r , (6) 
sin 23 IV + COs6l23*r, (7) 

which corresponds to multiplying the neutrino column vector in Eq. ([3]) with T23 from the left. Eq. ([3]) then becomes 




T13T12 




\t\ 2 tI 



V c ' 


.000, 




(8) 



1 Such findings are in agreement with what was found in Ref. |28l| . 

2 A remark on this aspect was made in [29| . 

3 The results we show in the present work do not include the difference in the mu and tau refractive indices which appear at one loop level 
due to different muon and tau lepton masses, V M r HI|], which is 1CT 5 V c . In fact, we have tested that the inclusion of this correction 
modifies very little the numerical results presented in Section III. 
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We define 

(Ex o o \ /v c o oy 

H = T 13 Ti 2 £ 2 T" l2 Tl + 0. (9) 
V E 3 J \ 0/_ 

The Hamiltonian depends on the CP-violating phase, (5. It is lengthy but straightforward to show that 

H(S) = 5+^(5 = 0)5, (10) 

where 

= ( 1 ) . (11) 




We are interested in solving the evolution equation corresponding to Eq. ([5]) : 

ih— = HU. (12) 
at 

It is important to recall that we need to solve this equation with the initial condition 

U(t = 0) = l. (13) 



Defining 

and using the relation in Eq. (fTU|) we get 



U = SU, (14) 



ih^=H(5 = 0)U , (15) 
at 



i.e., Uq provides the evolution when the CP-violating phase is set to zero. Using Eq. (fT3")) we see that the correct initial 
condition on U is Mo(t = 0) = S. However, Eq. (fH))) is nothing but the neutrino evolution equation with CP-violating 
phase set equal to zero. If we call the solution of this equation with the standard initial condition Uo(t = 0) = 1 to 
be Uo, we see that we should set U n = UoS, which yields 

U(6) = S^UqS. (16) 

Eq. (| 16[) illustrates how the effects of the CP-violating phase separate in describing the neutrino evolution. It is valid 
both in vacuum and in matter. This result is in agreement with (28l |. It is easy to verify that this result does not 
depend on the choice of the parametrization for the neutrino mixing matrix. 

Using Eq. (fT6|) it is possible to relate survival probabilities for the two cases with 6 = and 5^0. We define the 
amplitude for the process v x — > v y to be A xy when 5^0 and to be B xy when 6 = so that 

P(v x ^v y ,5^0) = \A xy \ 2 . (17) 

and 

P(u x ^u y ,S = 0) = \B xy \ 2 . (18) 

Using Eq. one can immediately see that the electron neutrino survival probability does not depend on the 

CP-violating phase. One can further write 

C23^4/ie — s 23^-re = C23-B Me — S2^B Te 
S2sA^ e + C2zA Te = e~ t5 [s2zB^ e + C2 3 B Te ] 

By solving these equations one gets 

A„ e = (c 2 23 + 4 3 e- u )B^ + c 23 S23(e^ - l)B Te , (19) 
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and 

Are = c 23 S23{e- lS - l)B^ e + (s 2 23 + c 2 23 e- iS )B re . (20) 

Clearly the individual amplitudes in Eqs. (fT5|) and (j2"0")) depend on the CP-violating phase. However, taking absolute 
value squares of Eqs. (fl"9|) and (|20|) . after some algebra one obtains that 

\A^ e \ 2 + \A Te \ 2 = \B^\ 2 + \B Te \ 2 , (21) 

or equivalently 

P{v„ ->v e ,S^0) + P(v T -> i/ e , 5 ^ 0) = P(i/„ -» zy e , 5 = 0) + P(^ r -» i/ e , 5 = 0). (22) 

Since P(v e — > i/ e ) does not depend on ^, Eq. (|22|) implies that if one starts with identical spectra with tau and mu 
neutrinos, one gets the same electron neutrino spectra no matter what the value of the CP-violating phase is (also 
see Eq l24]) . This result was first established in with a different derivation. A remark on this aspect is also made 
in 

Some differences in the muon/tau neutrino fluxes at emission can arise at the level of the Standard Model, from 
example from radiative corrections to the muon and tau neutrino cross sections |4lj . On the other hand, if physics 
beyond the Standard model operates during the infall and the shock-bounce stages of the supernova evolution, mu 
and tau neutrino fluxes can differ and induce CP-violating effects in the supernova environment. For example, generic 
neutrino-flavor changing interactions can give rise to significant net mu and tau lepton numbers [32l |. In particular, if 
there are flavor changing interactions involving charged leptons (e.g. a large scale conversion in the e~ — > fj,~ channel) 
one could also end up with significantly different mu and tau neutrino fluxes. In such cases one could have effects 
from the CP-violating phase on the electron (anti)neutrino fluxes as well. 

Note that our findings are at variance with those of Ref.[14(. In fact, there the authors conclude that even if mu and 
tau neutrino fluxes are different, CP-violation effects cannot be observed. Such a difference arises from the fact that 
different initial conditions are taken in our calculations compared to those used in Eq.(47) of fl4j. Indeed since the 
initial neutrino states should be those at the neutrinosphere, the neutrino conversion probability P ia should depend 
on the 8 phase (see section 3.4 of [HI). 



III. NEUTRINO FLUXES AND ELECTRON FRACTION IN SUPERNOVAE 



In this work we will discuss possible effects induced by the Dirac CP violating phase on two particular observables 
in the core-collapse supernova environment: the neutrinos fluxes <j) v and the electron fraction Y e . Note that the impact 
of the neutrino magnetic moment on such observables was studied in (33J. According to supernova simulations, the 
neutrino fluxes at the neutrinosphere are quite well described by Fermi-Dirac distributions [3J| or power-law spectra 
[35l | . Neutrino masses and mixings modify this simple pattern by mixing the spectra during neutrino evolution. Since 
muon and tau neutrinos only undergo neutral current interactions, they decouple deeper in the star. Electron (anti- 
neutrinos experience both charged and neutral current interactions, the anti-neutrino cross sections being weaker 
than for neutrinos and matter being neutron-rich. As a result a neutrino hierarchy of temperatures is expected, 
(E„ e ) < {E 9e ) < (E Ur ) with typical ranges of 10-13, 13-18 and 18-23 MeV respectively [35[. In order to show possible 
CP-violating effects on the Vi fluxes, we will use the ratio: 

where the neutrino fluxes are given by 

(t> Vi (S) = L Vi P{v t -> Vi ) + L V] P(v 3 -» + L Vk P(v k -> Vi) (24) 

with the luminosities 

L° E 2 
L Vi {r, E v ) - 4^2(^)3^)^) 1 + exp (eJt v _ v) ( 25 ) 

where i*2(*?) is the Fermi integral, Lq is the luminosity that we take as 6 x 10 51 erg/s as an example and and r is the 
distance from the proto-neutron star. We consider the Fermi-Dirac distribution as typical example. The quantities 
P{vi — > Vi/j) correspond to the survival/appearance neutrino probability during the evolution in matter. 
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The dominant reactions that control the neutron-to-proton ratio outside the hot proto-neutron star is the capture 
reactions on free nucleons 

i/ e +n#p + e~, (26) 

and 

i/ e + pv^n + e + . (27) 

We designate the rates of the forward and backward reactions in Eq. (f2"6"|) to be X Ue and A e - and the rates of the 
forward and backward reactions in Eq. (|2T[) to be A Pe and A e +. The electron fraction, Y e , is the net number of 
electrons (number of electrons minus the number of positrons) per baryon: 

Y e = (n e - - n e+ )/n B , (28) 

where n e - , n e +, and ng are number densities of electrons, positrons, and baryons, respectively. If no heavy nuclei 
are present we can write the rate of change of Y e as 

dY 1 

= A„ — (A p + X n )Y e + — (A p — X n )X a , (29) 

where we introduced the alpha fraction X a , the total proton loss rate X p — Ap e + A e - and the total neutron loss rate 
A n = X Ve + X e + . From Eq. (|29p one can write the equilibrium value of the electron fraction 

Ap + An 2 Ap + An 

As the alpha particle mass fraction increases more and more free nucleons get bound in alpha particles [3(|. This 
phenomenon, called alpha effect, pushes the electron fraction towards towards the value 0.5 (cf. Eq. (|50"|) ). Since it 
reduces available free neutrons, alpha effect is a big impediment to r-process nucleosynthesis [37| . At high tempera- 
tures, alpha particles are absent and the second term in Eq. (|30p can be omitted. Since electron and positron capture 
rates are very small, the electron fraction can be rewritten as 

e 1 + Ap/A„ ^ ^ 

with the capture rates on x = p, n given by 

A„,p = J a Uerii u eP (E v )cj)(E u )dE v (32) 
and a Vc n,v c p being the reaction cross sections for the corresponding processes Eqs. (|26H27]) . 



IV. POSSIBLE CP VIOLATION EFFECTS: NUMERICAL RESULTS 



It is the goal of this section to investigate numerically effects induced by the Dirac i5 phase : i) on the muon and 
tau neutrino fluxes when their fluxes at the neutrinosphere are supposed to be equal; ii) on the electron, muon, 
tau (anti)neutrino fluxes, when the muon and neutrino fluxes differ at the neutrinosphere. In fact, Eqs.(fTl |2"2"]) and 
([24]) show that in the latter case the electron (anti)neutrino fluxes become sensitive to the CP violating phase. We 
have performed calculations for several values of the phase. The effects discussed here are present for any value and 
maximal for 5 = 180°. For this reason most of the numerical results we show correspond to this value. 

We have calculated the neutrino evolution outside the supernova core using Eqs.([2]© and determined the neutrino 
fluxes Eqs. (|23H25|) and the electron fraction (|3H32[) . The numerical results we present are obtained with a supernova 
density profile having a 1/r 3 behavior (with the entropy per baryon, S — 70 in units of Boltzmann constant), 
that fits the numerical simulations shown in [24| . The neutrino fluxes at the neutrinosphere are taken as Fermi-Dirac 
distributions with typical temperatures of T„ e =3.17 MeV, T Pc =4.75 MeV and T Vx = 7.56 MeV (with v x = v ll .,v T , P M , D T ) 
(Figure QJ (the chemical potentials are assumed to be zero for simplicity). The oscillation parameters are fixed at the 
present best fit values [3j|, namely Amf 2 = 8 x 10" 5 eV 2 , sin 2 26» 12 = 0.83 and Am^ 3 = 3 x 10~ 3 eV 2 , sin 2 26» 2 3 = 1 
for the solar and atmospheric differences of the mass squares and mixings, respectively. For the third still unknown 
neutrino mixing angle 6*13, we take either the present upper limit sin 2 26'i3 = 0.19 at 90 % C.L. (L) or a very small 



6 



le+13- 




20 40 60 

Neutrino Energy (MeV) 



80 



FIG. 1: Neutrino fluxes at the neutrinosphere: the curves show the Fermi-Dirac distributions used for electron neutrinos 
with T„ e =3.2 MeV (solid), electron anti-neutrinos T Se =4.8 MeV (dashed) and for the other flavors T„ x = 7.6 MeV (with 
fx = Vfi, Vt, i/p, v T ) (dotted line). 
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Neutrino Energy (MeV) 

FIG. 2: (lower curves) and v T (upper curves) flux ratios for a CP violating phase 8 = 180° over 5 = 0° Eq. (|23[) . as a function 
of neutrino energy. Results at different distances from the neutron-star surface are shown, namely 250 km (dotted), 500 km 
(dashed), 750 km (dot-dashed) and 1000 km (solid line). The curves correspond to the normal hierarchy and sin 2 2#i3 = 0.19. 



value of sin 2 20i3 = 3x 10~ 4 (S) that might be attained at the future (third generation) long-baseline experiments 
Note that the value of 6*13 determines the adiabaticity of the first MSW resonance at high density [H, |4(| , while #12 
governs the second (adiabatic) one at low density. Since the sign of the atmospheric mixing is unknown, we consider 
both the normal (N) and inverted (I) hierarchy. In the former (latter) case (anti)neutrinos undergo the resonant 
conversion. We will denote results for the normal hierarchy and sin 2 2# 13 = 0.19 (N-H), inverted and sin 2 26* 13 = 0.19 
(I-H), normal hierarchy and sin 2 26»i 3 =3. 10~ 4 (N-S), inverted and sin 2 26»i 3 =3. 10~ 4 (I-S). 

Figures [2] and [3] show the v^, v T and v^, v T flux ratios Eq. ([23| for 8 = 180° over for 5 = 0°. One can see that large 
effects, up to 60 % are present for low neutrino energies in the anti-neutrino case; while smaller effects, of the order 
of a few percent, appear in the neutrino case. The effect of a non-zero delta over the v T fluxes as a function of 
neutrino energy is shown in Figure 2] at a distance of 1000 km. We see that an increase as large as a factor of 8 (4) 
can be seen at low energies in the (u T ) spectra. A similar behavior is found in the anti-neutrino case. 

In most of supernova simulations, the and v r luminosities are approximately equal, because these particles 
interact via neutral current only, at the low energies possible at supernova 4 . Since the v ei v e appearance probabilities 
are independent of 8 and as long as the and v T luminosities are taken to be equal, using Eqs. (|22[) and (|24[) one can 



4 Note however that, even at the level of the Standard Model, some differences can arise for example from loop corrections )4lH . 
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FIG. 3: Same as Fig[2]for (left) and v T (right) fluxes. 
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FIG. 4: Ratio of the (left) and v T (right) fluxes for S = 180° over S = 0° at a distance of 1000 km from the neutron-star 
surface. The curves correspond to N-L (solid), N-S (dashed), I-L (dot-dashed), I-S (dotted). 



show that the v e and v e fluxes are independent of the CP violating phase. Practically all the literature concerning 
the neutrino evolution in core-collapse supernovae ignore the Dirac phase, for simplicity. Our results justify this 
assumption if such calculations make the hypothesis that the (P M ) and v r (p T ) luminosities are equal and neglect 
the V^.r. 

On the other hand, the situation is different if the muon and tau neutrino fluxes are different at the neutrinosphcrc 
either because of the corrections within the Standard Model and/or because of physics beyond the Standard Model, 
such as flavor changing interactions [32| which might populate the , v T fluxes differently and differentiate their 
temperatures at decoupling. Our aim is to show the CP violating effects in this case. We have explored various 
differences between the ^ M ,^ T luminosities. We present here for example results corresponding to the v ^ , v T total 
luminosities Eq.flMD different by 10 %, e.g. L% t = 1.1 L° or T^= 8.06 MeV while T !/fj = 7.06 MeV. Figure [5] presents 
as an example the evolved v e ^v il neutrino fluxes, at 1000 km from the neutron star surface, when T„ ^ T Vt . The 
different curves show results for the two hierarchies and the two values of #13 . Similarly to the case where T„ = T„ T , 
while for the N-L case the first resonance is adiabatic and the electron neutrinos get a hotter spectrum, for all other 
cases the spectra keep very close to the Fermi-Dirac distributions (Figure [T]) . The situation is obviously reversed for 
the muon neutrino flux. Figures [6l [7] and [8] show the ratios of the v e and fluxes for a non-zero over a zero delta 
phase, as a function of neutrino energy, at different distances from the neutron star surface. One can see that effects 
up to a factor of 2-4 on the v e and up to 10 % on are present. A similar behavior is found for the v e and v T fluxes. 

The behavior of the flux ratios shown in Figure [7] is easy to understand. From Eqs. (fT^)) and (|20p one can write 



<M<*) = <^e (S = 0)+ sin 26*23 sin - (L Vt - L„J 



(33) 
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Neutrino Energy (MeV) Neutrino Energy (MeV) 



FIG. 5: Electron (left) and muon (right) neutrino fluxes Eq, (124[) at fOOO km from the neutron star surface, N-L (solid), N-S 
(dashed), I-L (dot-dashed), I-S (dotted). In the N-L case, the first resonance is adiabatic and the Fermi-Dirac distributions 
at the neutrinospere (Fig[TJ are completely swapped with u x . The situation is reversed for v^. These results are obtained by 
fixing T„ T larger than T„ by 1 MeV, as an example of the difference that could be induced by the presence of flavor-changing 
interactions in the neutrinosphere (see text). 
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FIG. 6: Ratios of the v e flux 8 — 180° over for 8 — 0° at 200 km from the neutron star surface, obtained by taking L% T = 1.1 L° 
(right) or T„ T = 8.06 MeV and T Vfl = 7.06 MeV (left) (see text). The curves correspond to N-L (solid), N-S (dashed"), I-L (dot- 
dashed), I-S (dotted). 



sin 2<? 23 sin - (| B Me | 2 - \B Te 



cos 26*23 sin - - % cos - {B^ e B* e ) + h.< 



Clearly the ratios calculated in these figures would be identity at the value of the energy where and v T spectra 
would cross (i.e., L„ T = L Ufl ). Away from this energy one expects an oscillatory behavior due to the additional terms 
in Eq. ([55)1 as the figure indicates. Note that even for 6 = the neutrino fluxes could also exhibit an oscillatory 
behavior. Concerning Fig. one can see that the effects due to S ^ and induced by taking different temperatures 
or luminosities are small, compared to the case with <5 ^ only (Figure H|) . 

Figure [5] shows results on the electron fraction Y e . Note that if 5 ^ there are no CP violation effects on Y e 
since this quantity depends on the electron neutrino and anti-neutrino fluxes only Eqs. (|31l32[) . Our results show that 
the effects due to S ^ are small (of the order of 0.1%) in all the studied cases with different muon and tau total 
luminosities and/or temperatures. 

Finally, we discuss the effects induced by the CP violating phase S on the supernova neutrino signal in a terrestrial 
observatory. Figure [TO] presents the expected number of events associated to electron anti-neutrino scattering on 
protons for different 5 values. This is calculated by convoluting the fluxes from Eq. (|24H25l) by the relevant anti- 
neutrino proton cross section [42j |. A water Cerenkov detector such as Super-Kamiokande (22.5 Ktons) is considered 
as an example. We assume 100 % efficiency. Note that the neutral current signal which is sensitive to all fluxes turns 
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FIG. 7: Same as Fig© but at 1000 km. 
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FIG. 8: Same as Fig[7]but for the i/ M flux ratios. 



out to be 5 independent as well, as can be shown by adding the three fluxes Eq. (j2~4")l . One cas see that S phase induces 
small modifications up to 5 % in the number of events, as a function of neutrino energy, and of the order of 2.10~ 4 on 
the total number of events. In fact, for a supernova at 10 kpc, we get for inverted hierarchy and large third neutrino 
mixing angle 7836.1 for 5 = 45°, 7837.0 for S = 135°, 7837.2 for S = 180°; while it is 7835.9 for 5 = 0°. These results 
are obtained with muon and tau neutrino fluxes having difference temperatures. Similar conclusion are drawn if we 
take different luminosities. For normal hierarchy and large 6*13, effects of the same order are found while for small #13 
and inverted/normal hierarchy the effecs become as small as 10~ 5 . 



V. CONCLUSIONS 



In this work we have analyzed possible effects induced by the CP violating Dirac phase in a dense environment such as 
the core-collapse supernovae. Our major result are that in matter: i) significant effects are found on the muon and tau 
neutrino fluxes for a non-zero CP violating phase; ii) important effects are also found on the electron (anti)ncutrino 
fluxes if the and v T neutrino fluxes differ at the neutrinosphcre. On the other hand the usual assumption of ignoring 
the CP violating phase made in the literature is justified if contributions from physics beyond the Standard Model is 
small and the and v T fluxes are equal at emission. We have calculated the events in an observatory on earth and 
shown that effects at the level of 5% are present on the number of events as a function of neutrino energy. 

Recent calculations have shown that the inclusion of neutrino-neutrino interaction introduces new features in the 
neutrino propagation in supernovae. A detailed study of the neutrino evolution with the CP violating phase, the 
neutrino-neutrino interaction as well as loop induced neutrino refractive indices will be the object of further work. 
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FIG. 9: Electron fraction for 5 = (left) and ratios of the electron fraction (right) for S = 180° compared to 8 = 0°, as a 
function of the distance from the neutron-star surface. The initial z/ M , v T fluxes have temperatures which differs by 1 MeV (see 
text). The results correspond to the normal hierarchy and sin 2 2#i3 = 0.19. 
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